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Abstract 

Starting from a second-quantized Hamiltonian of many-particle systems, we derive 
a set of time evolution equations for four-wave mixing (4WM) processes of coherent 
matter waves, which is analogous to those in optical 4WM except that the spatial 
variable is being replaced by time. Several interesting problems in 4WM such as the 
phase-matching condition, the effects of relative phase difference, and the conversion 
efficiency are then discussed in detail. We also show that the main features in recent 
4WM experiment [Deng et al, Nature 398, 218 (1999)] can be undersood within the 
present simplified model. 

PACS number(s): 03.75.Be, 03.75.Fi 



Introduction 



The experimental realizations of Bose-Einstein condensation (BEC) in trapped 
dilute alkali atomic vapours [1] and the atom "laser" [2] have made possible the 
experimental study [3] of four- wave mixing (4WM) of coherent matter- waves. 
As a typical effect of nonlinear atom optics [4], the possibility of 4WM with 
matter waves of Bose-Einstein condensate (BEC) can be deduced from the 
similarity between the Maxwell equations for lasing action in nonlinear optical 
media[5] and the Gross-Pitaevskii(GP) equation[6] widely used in studying the 
properties of macroscopic quantum systems such as BECs. 

In the past few years, 4WM in BECs and related problems have been studied 
intensively [3,7-16]. Early in 1995, Goldstein et al. had started to investigate 
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the problem of phase conjugation [7] in BECs, including the case of multi- 
component BECs[8]. A 4WM experiment using three colliding EEC wavepack- 
ets with different central momenta was proposed by Trippenbach et al. [9] . Re- 
cently, Deng et al. successfully carried out 4WM experiment in sodium conden- 
sates[3] by means of Bragg diffraction technology [17]. In addition, stimulated 
by the work of Ho and coworkers[10] and Ohmi and Machida[ll] on the Bose 
condensates with internal degrees of freedom, Law et al. studied quantum spin- 
mixing in spinor BECs [12] using an algebraic method developed in problems 
related to cavity QED[18], and Goldstein and Meystre[13] developed a quan- 
tum theory of atomic 4WM involving both the internal and spatial degrees 
of freedom. Rz§zeewski et al. investigated the fluctuations in the populations 
of atoms in 4WM wavepackets due to the quantum mechanical nature of the 
mean-field BEC wavefunctions[14]. Although the experimental results turn 
out to be consistent with the numerical simulations by Trippenbach et a/. [9] 
based on the GP equation, further theoretical studies are still needed in order 
to understand 4WM in BECs more clearly. Without invoking the undepleted 
pump approximation, Wu et a/. [15] calculated analytically the 4WM in the 
experimental condition of Ref. [3]. Based on the slowly- varying-envelope ap- 
proximation, Trippenbach et a/. [16] developed a three-dimensional quantum 
mechanical description for 4WM of wavepackets created from BEC. 

The main purpose of the present work is to study the ideal 4WM with co- 
herent matter waves from BEC in which only four plane waves with definite 
wavevectors and frequencies are involved. We organize the paper as follows. 
A set of time evolution equations for 4WM are derived in Sec. 2 starting from 
a Hamiltonian for many-particle systems in the second quantized form rather 
than the original GP equation. Section 3 gives a short discussion on the phase- 
matching condition for ideal 4WM. In Sec. 4, we present results demonstrating 
the effects of relative phase difference and the collapse and revival (C&R) 
behaviour of ideal 4WM. The problem of conversion efficiency and the main 
characters of recent 4WM experiment [3] are investigated in Sec. 5. Sec.6 gives 
a brief summary. 



2 Time evolution equations for 4WM 

To study 4WM with matter waves, we consider the dynamics of N interacting 
identical bosons. For simplicity, we omit the internal degrees of freedom. In 
standard experimental setups for studying 4WM, the confining potential is 
removed, and the Hamiltonian can be generally written in terms of bosonic 
creation and annihilation operators as [19] 

^ = IZ£fe4«fe+ ^^^LoL^fcaafei, (1) 

k k4+k3=ki+k2 
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where — Ti'^k'^ /2M is the kinetic energy of a particle in a single-particle 
state of momentum p = ^k, q = k4 — ki = k2 — ka denotes the momentum 
transfer in two-body scattering, and M denotes the atomic mass. Momentum 
conservation in the second term results from the calculation of the interaction 
matrix elements C/k4k3k2ki = C/<?(^k4+k3,ki+k2, where Uq = f (frU{r) exp(iq-r) is 
the Fourier transform of the interaction potential f/(r). The field operator can 
be expanded as ^(r, t) — J2k ^fc exp(ik • r) / ^/V with V denoting the volume of 
the system. 

Within the Heisenberg picture, the wavefunction is independent of time and 
the dynamical evolution of the system is given by that of the field operators. 
The annihilation operators satisfy equations of the form 

if^^^ = {£k4 + Uon)ak^ + Yl y43«fc2«fen (2) 

kl+k2-k3=k4 



where the subscripts ki,k2, and k^ can be any allowed values of wavevectors 
(or particle momenta), and the operator of the averaged number density of 
atomic gas is a fixed constant defined by n = (1/V) J2k olofe- 

Guided by results obtained through the experimental realization of BEC and 
the technique of Bragg diffraction [3, 17], we assume that the creation and an- 
nihilation operators a\ and of a certain mode with macroscopic population 
of atoms can be treated as a pair of conjugated complex numbers and ak, 
and those for the other modes vanished[ll,20,21]. Therefore the field opera- 
tor becomes a macroscopic wavefunction satisfying the GP equation, except 
that Uq may depend on q for scatterings with large momentum transfer [6, 22]. 
Introducing a new variable Ak = ak/VV and assuming for simplicity a q- 
independent Uq ~ Anaoh? /M with oq denoting the s-wave scatteing length, 
the slowly varying envelope function of the field, = A^ exp[i{ek — nUo)t/h] 
with n — Y^k^k^k/^ 1 satisfies a set of coupled nonlinear equations 



I 



dAk^ _2aoh ^ a* a a „iAu>t 



kl+k2-k3=k4 



E K^k^^k^e'^^' (3) 



with Au = (e^^ + £k3 ~£k2 The momentum conservation indicated in 

the summation is a consequence of the translational invariance of the system. 
In contrast, collision processes with Acu ^ are not forbidden in principle. 
As the system has a fixed number of atoms, the averaged number density is 
a constant, i.e., J2k l^fcP = "^^ We can then define the normalized field ampli- 
tudes as = ^k/V^j which satisfy the normalization condition J2k iCfcP = 1; 
to describe the relative distribution of atoms more conveniently. The time 
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evolution equations for these normalized amplitues take on the more compact 
form of 

= ^ E a.^k.^k.e^''^ (4) 



where the re-scaled time r = t/T is dimensionless with the characteristic time 
T = M/4naoh being inversely proportional to the number density, and the 
phase mismatch or 'energy-loss' parameter is defined as = TAlo. Before 
going into details, some important conclusions can be drawn on the common 
properties of 4WM with matter waves. The number density of atoms can affect 
the time evolution of 4WM only by modifying the characteristic time T, imply- 
ing that 4WM evolves faster in situations of higher number density of atoms. 
In actual experimental situations [3], the characteristic time is estimated to be 
about 3ms (see Sec. 5). 

The general form of time evolution equations for multi-wave mixing (MWM) is 
a rather complicated set of coupled nonlinear equations because the subscripts 
ki,k2,ks and in Eq.(4) can take on any allowed values of wavevectors (or 
particle momenta) . However it can be simplified for an ideal 4WM process to 

^ = -^e2*e3e4e-'"^ ^ = -tCiU,e-^''\ (5) 
^ = -^a6ae^"^ = -iCsUie^""'- (6) 



This set of equations is analogous to those in optical 4WM except that the 
spatial variable is replaced by the re-scaled time. In what follows, we will study 
4WM of matter waves based on the simplified model described by Eqs.(5) and 
(6). All numerical calculations are carried out using the 4th order Runge-Kutta 
method. 



3 Phase-matching condition for 4WM 



The recent 4WM experiment [3] shows that three colliding wavepackets with ar- 
bitrary central momenta cannot always create a fourth wavepacket, and a pos- 
sible 4WM process must satisfy the phase-matching condition, i. e. momentum- 
energy conservation. The Hamiltonian (Eq.(l)) describes a system satisfying 
the conservation of total energy, momentum and particle-number. Therefore 
the phase-matching condition is always understood as k4 + = ki + k2 
and -|- £fc3 = + ^ki ■ Note that from the dispersion relation of massive 
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particles Sk = {hky/2M, the phase mismatch parameter can be written as 
Q = (ks — ■ (k3 — ki)/87rnao. The phase-matching condition, i.e. fl — 0, 
for matter waves can be shown graphically as in Fig.l. 

It has also been emphasized, however, in the discussion after Eq.(l) that the 
scattering process of bosons does not require the conservation of energy. So 
the constraint of energy conservation deserves some more explanations. In 
fact, numerical calculations show that the parameter fl = TAcu affects 4WM 
processes, but only large enough phase mismatch can suppress the creation 
of fourth wave. Figure 2 shows the population A^4 of the created wave as a 
function of r for different values of Q. For small fl, N^i^r) is similar to that of 
an ideal 4WM process satisfying the phase-matching condition exactly. The 
conversion efficiency is lower for higher values of fl, and no 4WM signal can 
be observed when 1^21 > tt. 



4 Dynamics of Ideal 4WM 



The time evolution of an ideal 4WM system can be described by Eqs.(5) and 
(6). This kind of process satisfies not only particle-number conservation but 
also the conditions 

dr dr dr dr 



implying that the two pumping waves must provide atoms to amplify the 
probing and created waves at the same rate. When one of the two pumping 
waves is exhausted, the another one will take on the role of the probing wave. 
Figure 3 shows the calculated results for typical collapse and revival (C&R) 
behaviour of 4WM. Initially, the N atoms are distributed in three states cor- 
responding to the pumping wave 1, pumping wave 2 and probing wave 3. Due 
to two-body interaction and the bosonic enhancement effect, the atoms in the 
pumping waves begin to transfer into the states 3 and 4. After a short period 
of time, all the atoms are transferred out of state 2 and no more atoms are 
transferred into state 4. The second phrase of 4WM begins at around r ^ 5. 
Now waves 3 and 4 play the roles of the pumping waves, and wave 1 acts as a 
probing wave, wave 2 will become appreciable again until the pumping wave 
4 vanishes at about r 10.5. This completes one period of ideal 4WM. 

Another special dynamical character in 4WM processes is the effect of relative 
phase difference. It is not the phase of each individual wave that matters. 
Instead if we write = y^e**'', it is the combination of the phases a — 
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</'4 + V'a — </'2 — V'l ~ that affects 4WM in the way that 



dh 
dr 

da 



dr 



' dr 



dh 
"d^ 



1 



= 2^ sin a; 



-) cos a, 



(8) 
(9) 



where the amphfication factor is defined as ^ = \/hhhh- When one of the 
four waves vanishes, we can set the combined phase a to be an arbitrary value. 
Therefore the initial values of the pumping and probing waves do not affect 
the time evolution of atomic distributions. 



5 Conversion efficiency of 4WM 



All of the previous works[9,14,16] on the 4WM of matter waves have focussed 
on the initial rate of growth of the fourth wave because the time of wavepackets 
interaction in experiment [3] is very short (determined by the size and individ- 
ual momentum of wavepackets) . Initially the rate of growth of the fourth wave 
is proportional to the product of the intensities of pumping and probing waves. 
If the experimental setup can be changed in such a way that 4WM can pro- 
ceed for a much longer period, we need more precise knowledge about the 
conversion efficiency of 4WM. 

In the experiment [3] of 4WM with matter waves in BEC with N ^ 1.1 x 10^ 
sodium atoms (M = 3.848 x lQ~'^^g and = 28A), the pumping (1,2) and 
probing (3) waves are prepared by Bragg scattering[17]. It was estimated that 
the population numbers in the three states took on the initial values of 
6.8 X 10^ 6.9 X 10^ and 3.3 x 10^ respectively. Assuming that 4WM 

took place in a region of volume V — Li^ with L 0.1mm, the characteristic 
time will be T = ML^/ANaoh ^ 3ms. 

The four curves in Fig. 4 show typical results of the conversion efficiency 
N4/N as a function of the re-scaled time r obtained by Eqs.(5) and (6), hence 
assuming an ideal 4WM process with i7 = 0, for different sets of initial condi- 
tions. Curve (i) shows the situation in which = N2 ^ N^, which achieves 
a maximum efficiency of 50% if 4WM can proceed for a time long enough, 
e.g., t > lOT. Curve (ii) is calculated using a set of initial conditions similar 
to the experimental situation. Comparing with curve (i), it demonstrates a 
much greater initial rate of growth of the fourth wave and a, however, smaller 
saturation efficiency. Note that a conversion efficiency of 10% is achieved at 
t pa 1.9T pa 5.8ms, which is consistent with the results reported in Ref.[3]. The 
initial conditions for curves (iii) and (iv) are chosen so that both of them have 
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the same initial rate of growth of the fourth wave. Even so, the two curves 
demonstrate very different evolution processes as a function of r. Therefore, 
carefully prepared initial atomic distributions are crucial in achieving a higher 
initial rate of increase in the efficiency and at the same time a high enough 
saturation efficiency. 

While the initial conditions in Ref. [3] are close to those values in achieving 
a high efficiency, the observed value of 6 ~ 10% is still rather low. A higher 
conversion efficiency can then be achieved by increasing the interaction time 
and/or by shortening the characteristic time by increasing the number density 
of atoms. 



6 Summary 

In conclusion, a set of nonlinear equations for the time evolution of multi-wave 
mixing in BECs is dervied. The formalism leads to a physically transparent 
picture of 4WM with matter waves. The ideal 4WM model is successfully 
applied to analyse results of 4WM experiments. Factors affecting the value of 
the conversion efficiency are discussed. The GP equation in momentum space 
provides us with an efficient way to study the wavepacket effects of 4WM 
because the four wavepackets are limited to four spots of the momentum space 
even though they may spread widely in real space. 
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FIGURE CAPTIONS 



Fig.l The momentum sphere of phase- matching condition for 4WM. The 
center is located at (ki -|-k2)/2, and the radius of the sphere is defined as the 
magnitude of the vector (ki — k2)/2. If the momentum of the probing wave 
points exactly at a point on the surface of the sphere, the largest conversion 
efficiency of 4WM is achieved. 

Fig. 2 The conversion coefficiency as a function of interaction time for 
five different values of the phase- mismatch paramter Q — TAcv. The time is 
given in units of the characteristic time as defined in the text. The five curves 
from top to bottom correspond to Q/ir = 0, 0.1, 0.2, 0.4, 1, respectively. 
The initial conditions for the pumping and probing waves are chosen to be 
e?, ^2, ^1 = 0.68, 0.58, 0.45. 

Fig. 3 Typical time evolution curves of the populations of the matter 
waves involved in a 4WM process. The initial conditions are chosen to be 

e?, ^2°, = 0.68, 0.58, 0.45. 

Fig. 4 Typical results of the conversion efficiency as a function of the inter- 
action time for different initial conditions. The initial conditions for curve (i) 
are C?, = 0.707, 0.703, 0.243, which lead to a maximum efficiency of 

50% but with a low initial increasing rate. Curve (ii) has initial conditions 
similar to those in experimental situation: = 0.637, 0.632, 0.441. 

Curves (iii) and (iv) show that in order to achieve an appreciable efficiency, 
the probing wave can be weak but the pumping waves must be strong. 
The initial conditions are: (in) = 0.804, 0.553, 0.217 and (iv) 

e?, ^2°, ^3° = 0.804, 0.217, 0.553. 
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